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SUMMARY 


/ 

This  report  is  made  up  of  two  parts.  In  the  first  part  we  present  additional 
results  on  the  performance  of  random  access  protocols  for  mobile  packet  radio 
networks.  The  effect  of  three  different  types  of  diversity  reception  in  Nakagami 
fading  environment  on  the  channel  throughput  and  the  average  packet  delay  of 
nonpersistent  carrier  sense  multiple  access  (NPCSMA)  protocol  is  considered. 
Expressions  for  the  probability  of  packet  error  with  diversity  receptions  are  newly 
obtained  for  both  independent  and  correlated  diversity  branches.  A  noncoherent 
frequency  shift  keying  modulation  scheme  is  assumed. 

In  the  second  part  a  comprehensive  study  of  the  problem  of  synchronization 
over  fading  dispersive  channels  is  presented.  Synchronization  is  a  fundamental 
problem  in  digital  communications  as  used  in  mobile  radio.  In  the  first  chapter, 
we  consider  a  simple  binary  detection  problem.  The  effect  of  receiver  mismat  ch 
is  investigated.  A  closed  form  expression  for  the  probability  of  false  alarm  and 
probability  of  detection  are  given,  also  a  set  of  curves  are  provided  to  demon¬ 
strate  the  amount  of  degradation  for  under-spread  channels  with  some  special 
scattering. 

Chapter'll  deals  with  the  performance  of  serial  synchronizer  over  fading 
dispersive  channels.  The  performance  indices  are  the  mean  and  variance  of  syn¬ 
chronization  time.  Upper  and  lower  bounds  on  the  mean  and  variance  of  syn¬ 
chronization  time  for  a  very  general  serial  search  system  are  derived  and 
evaluated.  It  has  been  found  that  the  results  are  highly  dependent  on  the  scatter¬ 
ing  function  of  the  channel,  the  degree  of  the  spread,  number  of  the  cells  used  in 
the  search  procedure,  and  the  time  required  to  reject  an  incorrect  cell  when  a 
false  alarm  occurs.  An  optimum  closed  loop  structure  for  symbol  synchronizer  is 
derived  in  chapter  III.  The  optimum  synchronizer  is  similar  to  those  already 
known  symbol  synchronizers  that  are  being  used  over  additive  white  Gaussian 
noise  (AWGN)  channels,  except  that  the  control  signal  is  completely  random; 
even  in  the  absence  of  the  AWGN.  It  is  shown  that  the  synchronizer  structure  is 
highly  dependent  on  the  scattering  function  of  the  channel.  For  a  Gaussian 
shaped  scattering  function,  a  simple  closed  loop  structure  is  obtained  and  expres¬ 
sions  for  tracking  error  statistics  are  derived. 

In  chapter  IV  the  method  of  stochastic  approximation  is  applied  to  the  prob¬ 
lem  of  synchronization  over  fading  dispersive  channels.  A  recursive  estimation 
procedure  is  developed  to  estimate  the  two  parameters  needed  by  the  synchron¬ 
izer.  The  result  is  a  closed  loop  in  which  the  component  of  the  error  signal  is 
proportional  to  the  expected  value  of  the  derivative  of  the  likelihood  function 
with  respect  to  the  appropriate  parameter. 
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PERFORMANCE  EVALUATION  OF  DIVERSITY 
IN  RANDOM  ACCESS  PROTOCOLS 


I.  Introduction 

The  effect  of  Rayleigh  fading  on  the  throughput  of  different  random  access 
protocols  such  as  Aloha  and  NPCSMA  have  been  introduced  [1,2].  However,  the 
Nakagami  distribution  [3]  covers  a  wider  class  of  fading  than  any  other  model. 
Experimental  as  well  as  theoretical  [4,5]  studies  showed  that  it  is  the  best  fit  for 
data  obtained  from  a  model  of  urban  radio  multipath  channels.  The  effect  of 
diversity  reception  on  the  probability  of  bit  error  with  Nakagami  fading  was  first 
considered  in  [6].  The  analysis  was  confined  to  CPSK  which  is  not  suitable  for 
mobile  communications,  beside  that  these  results  cannot  apply  to  packet  switch¬ 
ing  where  we  have  packets  of  N  bits.  Some  other  results  were  obtained  [7,8]  for 
NCFSK  with  MRC  and  selection  combining  diversity  techniques.  Again  these 
results  were  obtained  for  the  probability  of  bit  error  as  well  as  the  third  combin¬ 
ing  technique,  namely  equal  gain  combining,  was  not  considered.  The  effect  of 
Nakagami  fading  on  the  NPCSMA  throughput  was  first  considered  in  [9]  for  suc¬ 
cessive  correlated  packet  retransmissions.  However,  the  analysis  of  the  diversity 
reception  on  the  probability  of  packet  error  and  consequently  the  channel 
throughput  represents  an  uninvestigated  case.  In  section  II,  the  probability  of 
packet  error  with  MRC  receiver  is  obtained  when  the  diversity  branches  are 
independent  with  either  identical  parameters  on  each  branch  or  not.  The  effect 
of  correlation  among  two  diversity  branches  is  also  considered. 

Similar  results  are  presented  in  Section  III  for  selection  combining.  Section 
IV  includes  the  results  for  equal  gain  combining.  Modified  expressions  for  the 
throughput  and  the  delay  are  obtained  in  Section  V  for  NPCSMA  protocol  along 
with  conclusions. 


II.  Probability  of  block  error  of  MRC 
diversity  system: 

The  following  assumptions  are  made  to  obtain  the  results: 

a)  NCFSK  modulation  mode  is  considered. 

b)  The  fading  on  each  branch  is  nonselective  flat  type  of  fading  that  is 
modeled  using  Nakagami  distribution. 

c)  The  noise  is  additive  white  Gaussian  with  spectral  height  2N0,  and  is  the 
same  at  all  branches.  One  of  the  known  diversity  techniques  [10]  is  the  MRC 
which  allows  the  signal  to  noise  ratios  to  be  added  noncoherently  along  the  M 
branches.  In  other  words 
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M 

E  Tk 

k-1 


(1) 


The  fading  envelope  S(t)  at  the  1th  branch  has  a  Nakagami-m  distribution  as  fol¬ 
lows: 


f(S)  = 


s^‘  "p(-ws) 


(2) 


where  =  <Sj2> 


.  m,  = 


Oj2  m>i 

<(si2-nj)2>  ,mi  2 


(3) 


The  signal  to  noise  ratio  at  the  ’1th  branch  will  be 


%  = 


2N„ 


(4) 


The  probability  density  function  (p.d.f.)  for  %  will  be 
(kif 


r(m.) 


(5) 


where 


2N-m,  m.  _ 

kj  ==  — — —  =  — -  ,  t j  =  Average  SNR 

“i  li 


(6) 


This  can  be  recognized  as  a  Gamma  distribution.  The  probability  of  packet  error 
is  defined  as  follows 


00 


Pp.  -  /  [1  -  (1  -  Pbfl  Hi)  <*1 


(?) 


where  N  is  the  number  of  bits  in  each  packet.  Pb  is  the  probability  of  bit  error 
for  a  certain  modulation  scheme.  For  NCFSK  Pb  is  given  [12]  by: 


Pb  =  y  . 1i  >  0 


(8) 


Case  (A)  :  Independent  Diversity  Branches 

For  identical  parameters  m,  0  on  all  branches,  we  have 
N  00 

Pp,  =  1  -  E  (q)  (4p  /  «~,/J  Hi)  i-t 


A  n  2 
n— 0  *  o 


(0) 


The  sum  of  i.i.d.  Gamma  distribution  with  parameters  m,k  is  also  a  Gamma  dis¬ 
tribution  [13].  i.e. 

t(l)  =  '7in'"1  exp(-k'y)  (10) 


where 
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mt  =  Mm 

Substituting  for  (10)  in  (0)  and  performing  the  integration  we  get 


' n'^  o'  n 

— «  2  (1  +  -i- 

'  2k 


Equation  (12)  can  be  easily  evaluated  when  N  is  not  large.  This  is  a  reasonable 
assumption  (14],  since  a  large  N  would  cause  an  unacceptable  probability  of  error. 
If  the  parameters  differ  in  each  branch  we  have 


Pp.  -  1  -  E  (, )  (4)  L-T-KM).  _  „ 

n=0  *  °  a 


where  L.T.  is  the  Laplace  transform  operator.  The  characteristic  (c/c)  function 
for  %  is  given  by  (13) 

“  <  K—r’"'  <“> 


Then  the  c/c  function  for  -7  will  take  the  form: 
M  K; 

_  TT  (  1 


-  £  <k4t> 


Consequently,  the  L.T.  for  f^)  will  be 

F(s)  =  &  >" 


Substituting  for  (16)  in  (13)  one  gets: 


N  vr  1  M  K; 

n  /  1  \n  n  f  1 


Pp.  =  1  -  E  <")  (~r  n  ( 

_ _ n  1=1 


K.  +  f 


For  arbitrary  branch  parameters  mk,  flk  but  equal  (77-)  for  all  branches  we  get 


L/7)  ~  Gamma  (k,  mk) 
k-l 


And  consequently 


,NW  1 


u 

E  “k 


Pp.  =  i-EQ(-f)”n(¥fS7jr 


Case  B:  Correlated  Branches  (M  =  2) 

In  this  case  we  assume  that  the  fading  envelope  has  the  same  parameters  in 
both  branches.  The  envelope  in  this  case  will  take  the  form  [3]: 


.w.v.v; 


-4- 


2\/jt  ,  m  .jm _ S  ,  -mS2 


f(S)  =  iBr  (^) 


r(m)  n  (i -p)m 


exp( 


n(W) 


2mvp  ®  -  y  n(l-p) 


(20) 


where  Ix(‘)  is  the  modified  Bessel  function  of  the  first  kind  of  order  X.  Substitut¬ 
ing  for  equation  (4)  in  (20)  one  gets:  < 

1 '  r(m)(i-p)ra  1  n 1  p(  n(i-p) 1 


,  n(i^X*N0-») ,» -  i  , 

(  ft.#  *- 


»-l(  fl(l-p) 


2mN0\/p/7 


) 


(21) 


Hence  the  probability  of  packet  error  will  take  the  form: 


N 


Ppe  =  i-c  e  0(-ir 

n=0  * 


oo 


f  ,  n^  x  ,  2N0C17 
/  exp(-  -£-)  exp( -^— ) 


where 


2N  Tf  m  —  — - 

2  *  1  (2No  Cl  7)  *y 

2^i  m - T 


C  = 


2N0>A 


m  \2m 


(— ) 

r(m)  (1  n 


(22) 


(23) 


,c,= 


m 


(x  -  />)( 1 

Rearranging  (22),  one  obtains: 


(24) 


—  N0  m  -  — *  N  xt  i 

Pp.  =  1  -  K(-^)  2  E  ©  (-  I)" 

U1  n-0  1 


\ 

i 


I 


'■'I-Vy-v 


■  -  •  *  * 
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oo 


/  (iY~1  Ix(b'>)  exp(-a7)d7 


(25) 


where  x  =  m-  i-,y  =  m  +  i-, 


n  2N0C.  , 

a=f +  ~^r*b  =  2N”Ci 


(26) 


The  integration  in  equation  (25)  is  found  to  be  [15] 


F(2m)  •  (2N0C,) 


-i 


r(m  +  i)2*-M|.  +  ^-r 


,  2F,(m  +  ~  ,  m  ,  m  +  -J-  f  (- 


2N0C! 


f  +  2N0C,/v^ 


(27) 


where  ,  is  the  hypergeometric  function. 

The  conditions  to  be  satisfied  for  the  above  integration  to  exist  are 


(1)  x  +  y  >  0 

(2)  |a|  >  |  b  | 

(3)  Re  |  a  |  >  Re  |  b  | 


(28) 


It  is  clear  that  equation  (28)  is  satisfied  for  every  value  of  n  and  m.  Using  the 
relations  from  [16] 


one  gets  for  P,*  : 


i  jF i(a>b,c,z)  —  t  i(b,a,c,z) 

,  2Fi(b,a,a,z)  =  (“  f3 

,  r(2m)  =  r(m)  r(m  +  j) 


Pr.“l-( 


2N0m 

(Wi-p 


\2m 


N  m 
£  (o>(- 

n— 0 


(20) 

(30) 

(31) 


1 


2N0C, 

~7T 


f  -  (2N0C,)r 


(32) 


* 
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III.  Probability  of  packet  error 

for  selection  combining: 

For  two  diversity  branches  we  have: 

Case  A:  Independent  fading  on  the  diversity  branches. 

The  resulting  signal  to  noise  ratio  can  be  described  by: 

1  =  Max  (7„  Tf2) 

Then,  for  arbitrary  parameters,  we  get: 

ffr)  -  f(li)  I ,  F^s)  +  f(lj)|11_1  F^,) 


Then 


f(i)  -  k,”1  «*p(-k,i) 


T  IS  m* 
r  m.-l 


r(m2) 


72  *  exp(-K272)  <b 2 


+  K"' 


2  r(m2) 

K,~ 


1  ®*-i 


7mrl  exp(-K27) 


‘  1'a'1  exP(~K^i) 


Consequently,  the  probability  of  packet  error  will  be: 

H,®'  K2®*  N  N  i 

U7 


if  f  *XP(~  ~  Kn) 
0  0  ^ 


■  exp(-K272)  *r2m' '  d-r2 

+  /  /  «p(-^  -  K2i)  ■T’-' 

o  o  ^ 

•  expC-Kjif!)  7®1"1  d  7\  d7}) 


In  the  above  integrals  let  x  =  7J7  ,  i=l,2,  Then  we  get: 


K,®'  K2“* 
r(m,)  r(m2) 


{/  x"*-1 


O 


/  exp  -  t(|-  +  Kt  +  K2x) 


-7- 


i  +  njf-l 


1  oo 

f  .“rl  f 


d7dx  +  /  X®*"1  / 


exp  '  l(f  +  Kj  +  Kjx)  V"' +  “^’d-rdx}] 


Or  using  the  definition  of  the  hypergeometric  function  and  making  some  manipu¬ 
lations,  we  arrive  at: 

N,  N 

B(mi,m2)  £,V  2ryn> 


r  1  1 

*  (—  *  — - 


ml  (  “.  +  K2)“‘ +  “* 
2 


JFiif&i  +  m2,  mi,  mx  +  1, 


K,+  f 
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aF^mi  +  m2,  m2,  m2  +  1, - )} 
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Case  (B):  Correlated  Branches : 

The  joint  pdf  of  two  SNR’s  7l,  7a  whose  envelopes  are  jointly  random  vari¬ 
ables  is  given  by  [9j: 

m-I 

f(7i>  Ha)  =  K'(7i72)  2 

•  exp  -  -p-  (7l  +  72)  In»-i(K“v'V72)  (39) 

1  -P 


where 


K‘  = 


K”*1  ,  K„  =  2^k 

“"I  ’  W 


r(m)(i-pX/»)  - 

Because  of  the  symmetry  of  7l  7a  in  (39),  one  can  write  for  f(7)  : 

t{i)  =  2  /  f(7,7i)d7i 


—  N  Kt  1 

Ppt  =  1  -  2  £  (£)(“-)“ 
n-0  z 


Hence 


•  1 1  «P(~7r)  Kwi)  *7  (42) 

o  o  ^ 

Lei  t®  =  7j  in  (42)  and  change  the  order  of  integrations  we  obtain  using  (15): 
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(43) 

It  is  obvious  that  the  conditions  for  the  integration  to  exist  are  satisfied 

V0<x<l,0<  |i|  <1. 

Using  a  similar  approach  to  case  A  and  after  some  simple  manipulations  one 
arrives  at: 

N 

^n' '  2 

/ _ »”~‘|ax  +  bl _ ^ 


where 
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L  =  2  r(2m) —  (  JiLyn 

r(m)T(m)  '  1-p ' 


a  == 
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The  above  integration  can  be  easily  evaluated  Vn  numerically. 


(44) 


(45) 

(46) 


IV.  Probability  of  Packet  (block)  Error 
for  Equal  Gain  Combining 

For  M  diversity  branches  with  identical  parameters,  the  fading  envelope  can 
be  described  by 

S(t)  =  £  Sk(t)  (47) 

k-l 


Case  A  Independent  diversity  branches: 

From  |3],  the  p.d.f.  of  S(t)  can  be  approximated  by: 


(48) 


where 
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-  MJfi(l  -  -21) 
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The  probability  of  packet  error  will  take  the  form: 
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=  Tf  and  integrate,  we  get 


Pp. = >  -  e  („)  (-  b" 
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,  Nx  =  MN0 


Case  B:  Correlated  diversity  branches  (M  =  £) 


We  assume  that  the  two  branches  have  identical  parameters  (m,ft).  The 
p.d.f.  of  the  combined  output  will  be 

s 

f(S)  =  /  f(S„  S-S,)  dS,  (53) 


Hence, 
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If  we  let  R|  =  Rx,  we  get 
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where  L  and  a  are  constants, 
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E±  K  n ’  ’  0(1  -p) 

__  r(m)(l-p)/»  2 

Using  (15],  and  after  satisfying  and  regularity  condition^  wTget 
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\£  Modified  Expressions  for  the  Channel 
Throughput  for  NPCSMA  Scheme 

The  channel  throughput  and  the  average  packet  delay  for  NPCSMA  proto¬ 
col  are  given  by  (17) 


S  =  G  exp(-aG)  \ 

G(1  +  2a)  +  exp(-aG)  I  ^ 

D  St  a  +  1  +  2a  +  (-|-  -  1)(1  +  2a  +  6+  a)  I 

The  modified  expressions  for  S  and  D,  if  the  packet  and. acknowledgement  chan¬ 
nels  are  fading  and  identical,  will  be 

„  '  ]  m 

Dn,Sa+l  +  2a+(f--l)(l  +  2a  +  Ha) 


For  the  last  3  diversity  schemes,  we  can  always  compute  the  expressions  for  Sm 
and  Dm  to  investigate  the  effect  of  different  combining  techniques  on  the  channel 

throughput.  _ 

Numerical  Results 

The  improvement  in  the  throughput  and  the  delay  of  a  channel  with  Naka- 
gami  fading  using  MRC  diversity  technique  is  displayed  in  figures  1  through  8.  It 
is  obvious  that  the  performance  approaches  that  of  a  nonfading  channel  with  4 
diversity  branches.  However  the  approach  is  quicker  when  the  branches  have 
different  parameters.  Correlation  among  two  diversity  branches  is  shown  to 
greatly  affect  the  performance.  Consequently,  serious  attention  must  be  directed 
to  assure  independence  between  branches.  It  is  also  noticed  that  the  limiting  case 
where  p  — *■  1  is  achieved  for  MRC  at  a  SNR  3  dB  less  than  that  needed  for  a  sin¬ 
gle  channel.  This  is  expected  from  equation  (1).  Similar  results  are  presented  for 
both  SGC  and  EGC.  For  SGC  the  performance  degrades  to  reach  that  of  a  sin¬ 
gle  channel  at  p  — ►  1.  For  SGC  the  limiting  case  as  p  — ►  1  occurs  at  an  average 
SNR  1.5  dB  less  than  that  needed  for  a  single  channel.  It  is  also  to  be  noticed 
that  for  SGC  a  high  SNR  branch  along  with  the  worst  fading  figure  (m),  the  per¬ 
formance  is  better  than  that  of  two  identical  branches  with  relatively  good  fading 
figure.  This  implies  that  the  average  SNR  affects  the  performance  in  a  more 
effective  way  than  the  fading  figure  (which  describes  the  severity  of  fading)  does. 
The  MRC  technique  shows  a  better  performance  followed  by  EGC  and  then  the 
SGC.  However,  the  improvement  in  performance  of  MRC  over  EGC  is  not 
attractive  to  favor  MRC  over  EGC  keeping  in  mind  that  EGC  is  easier  to  imple¬ 
ment  than  MRC.  The  results  obtained  here  contain  the  case  of  Rayleigh  fading 
(at  m  =  1)  as  well  as  the  case  of  single  bit  transmission  (at  N  =  1).  _________ 

VI.  Conclusions: 

The  probability  of  packet  error  for  three  different  diversity  combining  tech¬ 
niques  is  evaluated  for  independent  fading  envelopes  over  each  branch.  The 
effect  of  correlation  among  two  different  branches  was  also  considered.  The  fading 
was  modeled  according  to  the  Nakagami-m  distribution.  The  results  were  used 
to  obtain  modified  expressions  for  the  channel  throughput  and  the  average  packet 
delay  for  NPCSMA  protocol. 
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PART  n 

SYNCHRONIZATION  OVER  FADING  DISPERSIVE  CHANNELS 

CHAPTER  I 

PERFORMANCE  DEGRADATION  DUE  TO  RECEIVER  MISMATCH 
IN  COMMUNICATION  OVER  FADING  DISPERSIVE  CHANNELS  f 

The  engineering  importance  of  fading  dispersive  channels  has  increased 
markedly  in  the  recent  years.  Optimum  detection  over  such  channels  requires  a 
perfect  knowledge  of  the  channels  scattering  functions  which  completely  charac¬ 
terize  the  behavior  of  fading  dispersive  channels.  Lack  of  such  knowledge  may 
result  in  performance  degradation.  The  objective  of  this  study  is  to  analyze  the 
performance  degradation  under  mismatch  conditions.  Two  types  of  receiver 
mismatch  will  be  considered.  First,  the  receiver  is  assumed  to  have  knowledge  of 
the  shape  of  the  scattering  function  but  the  mean  time  delay  and  the  mean  fre¬ 
quency  shift  are  unknown.  Secondly,  we  will  investigate  the  performance  degra¬ 
dation  when  the  shape  of  the  scattering  function  is  not  known  to  the  receiver.  A 
closed  form  expression  for  the  probability  of  false  -  larm  and  probability  of  detec¬ 
tion  are  given,  also  a  set  of  curves  are  provided  to  demonstrate  the  amount  of 
degradation  for  under-spread  channels  with  some  special  scattering  functions. 

1.1  INTRODUCTION 

During  recent  years  there  has  been  an  increasing  amount  of  attention 
given  to  the  study  of  fading  dispersive  channels,  channels  that  exhibit  both  fad¬ 
ing  and  dispersion.  The  determination  of  optimum  modulation  and  demodulation 
techniques  and  the  evaluation  of  the  efficacy  of  optimum  and  suboptimum 
receivers  for  such  channels  are  of  major  importance  to  the  system  designer. 

Fading  dispersive  channels  are  usually  best  described  as  random,  linear, 
time-invariant  filters  [1  ].  The  characterization  of  time-variant,  linear  filters  in 
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terms  of  system  i  unctions  received  its  first  general  analytical  treatment  by 
Zadah  [2],  who  introduced  the  time-variant  transfer  function  and  the  bi-frequency 
function,  as  frequency  domain  methods  of  characterizing  time-variant  linear 
filters.  There  is  another  common  approach  to  describe  the  channel,  an  approach 
which  involves  the  notion  of  scatterers.  This  approach  leads  to  a  physical  picture 
of  the  channel  as  a  continuum  of  moving  scatterers.  In  such  model,  propagation 
is  established  by  a  single  scattering  from  a  large  number  of  independent  scatter¬ 
ers. 

In  many  applications,  it  is  reasonable  to  suppose  that  the  impulse 
response  of  the  filter  is  a  sample  function  of  a  Gaussian  random  process.  Using 
such  supposition,  the  specification  of  the  channel  reduces  to  the  specification  of 
the  mean  and  correlation  functions  either  of  the  channel’s  random  impulse 
response  (first  approach)  or  of  the  received  process  conditioned  upon  the 
transmitted  waveform  (second  approach).  Many  channels  are  adequately  modeled 
by  taking  the  mean  to  be  zero  and  the  correlation  function  to  be  a  special  form, 
which  is  determined  by  the  scattering  function  of  the  channel  as  will  be  seen  in 
the  next  section. 

In  section  1.2  mathematical  model  of  the  channel  is  presented,  and  the 
statistics  of  the  output  of  the  channel  are  given  in  terms  scattering  function  of 
the  channel  and  the  complex  envelope  of  the  transmitted  waveform. 

In  section  1.3  we  present  the  receiver’s  structure  and  discuss  a  way  of 
obtaining  an  approximation  to  the  probability  of  detection  and  the  probability  of 
false  alarm.  The  discussion  in  this  section  is  limited  to  special  category  of  prob¬ 
lems;  namely,  the  low-energy-coherence  problems. 

In  section  1.4  we  investigate  the  degradation  in  performance  due  to  the 
mismatched  receiver,  by  that  it  is  meant  that  the  receiver  has  riot  enough 
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information  regarding  the  scattering  function  of  the  channel.  The  two  cases  of 
interest  are 

(1)  The  shape  of  the  scattering  is  known,  except  for  some  unknown  parame¬ 
ters. 

(2)  The  scattering  function  of  the  channel  is  completely  unknown  to  the 
receiver. 

1.2  CHANNEL  MODEL 

The  channel  to  be  considered  is  depicted  in  Figure  1.1.  This  channel  is 
completely  described  by  the  scattering  function  cr(r,f).  The  most  important 
parameters  of  o(r,f)  are  B,  the  frequency  interval  in  f  outside  of  which  er(r,f)  is 
essentially  zero,  and  L,  the  time  interval  in  r  outside  of  which  cr(r,f)  is  effectively 
zero.  The  quantity  B  is  called  the  Doppler  spread,  and  represents  the  average 
amount  that  an  input  waveform  will  be  spread  in  frequency,  while  L  is  known  as 
the  multipath  spread,  and  represents  the  average  amount  in  time  by  which  an 
input  signal  will  be  spread.  The  product  S=BL  is  called  the  spread  of  the  chan¬ 
nel,  channels  for  which  S  <  1  are  called  under-spread  channels,  while  channels 
with  S  >  1  are  called  over-spread  channels.  Characterization  and  classification  of 
fading  dispersive  channels  are  given  in  [3]. 

Throughout  the  discussion,  we  denote  the  channel  input  and  output 
waveforms  by  s(t)  and  y(t),  respectively.  These  waveforms  will  usually  be 
represented  by  their  complex  envelopes,  thus 


s(t)=Re|u(t)exp  [ja;0t  J  j, 

(1.1-a) 

y(t)=Re|v(t)exp  [ju;0t  j  j, 

(1.1-b) 
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■where  Re  {.}  denotes  ”  the  real  part”  of  the  indicated  quantity,  and  denotes 
the  nominal  carrier  frequency  in  radians  per  second. 

In  many  applications  it  is  reasonable  to  assume  that  the  received  signal  is 
-  conditionally  Gaussian.  By  conditionally  Gaussian,  we  mean  that  given  the 
transmitted  waveform,  the  received  waveform  is  Gaussian  random  process.  If  one 
accepts  this  supposition,  the  description  of  the  channel  reduces  to  the 
specification  of  the  mean  and  correlation  function  of  the  received  process  condi¬ 
tioned  upon  the  transmitted  waveform.  It  can  be  shown  [3]  that 


Ry(t,r)=Re|R(t,r)exp  [ju>0(t-r)  j  j, 


where 


R(t,r)=  J  J  o(r,f)u(t-r)u*(r-r)exp  |j27rf(t-r)  jdrdf. 


The  function  R(t,r)  defined  by  (1.3)  is  the  complex  correlation  function  of  the 
complex  envelope  of  y(t),  that  is 


t,r)=E|v(t)v*(r)|. 


Insight  into  the  scattering  function  and  methods  of  measuring  it  are  discussed  in 
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1.3  RECEIVER  STRUCTURE 


Optimum  receiver  structure  for  digital  signaling  over  fading  dispersive 


channels  have  been  determined  [5],  and  the  performance  of  binary  signaling  sys¬ 


tems  has  been  thoroughly  investigated.  In  this  section  a  brief  discussion  of  the 


optimum  receiver  structure  for  simple  binary  detection  problem  is  presented.  The 


received  waveforms  under  the  two  hypotheses  are 


l.»V 


Hi: 

*(t>»v(t)+w(t), 

T,<t<Tf, 

Ho: 

|(t)=sw(t), 

T,<t<T„ 

(1.5) 

We  assume  that  that  w(t)  is  a  complex  white,  zero-mean  Gaussian  process  with 
spectral  height  N0  and  v(t)  is  zero-mean  complex  Gaussian  random  process  with 
covariance  function  R(t,r). 

The  approach  to  design  the  optimum  receiver  is  analogous  to  the  approach 
used  in  the  deterministic  signal  case  [6].  It  can  be  shown  that  the  optimum 
receiver  compares  the  likelihood  ratio,  l,  with  a  threshold  tj, 

I  <  Vi  (1.6) 

Hi 

where 

j  T' 

z*(t)h(t,r)z(r)dtdr,  (1.7) 

T| 

and  h(t,r)  satisfies  the  integral  equation 

T, 

N0h(t,r)+Jh(t,u)R(u,r)du=R(t,r),  T,<t,r<T(.  (1.8) 

T, 

The  solution  to  (1.8)  can  also  be  written  in  terms  of  eigenfunctions  and  eigen¬ 
values  of  the  complex  correlation  function  R(t,r)  (8), 

hM-E  x,+No/2  *.(t)*'(r).  <>■») 

One  possible  realization  to  generate  l  is  the  filter-correlator  receiver  see  [5]. 
Even  though  the  procedure  is  well  defined,  the  actual  implementation  is  difficult, 

mainly  because  of  the  dependence  of  h(t,r)  on  X,  and  ^,(t);  i= 1 ,2 . There  are 

several  categories  of  processes  for  which  one  can  obtain  a  reasonably  nearly 
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optimum  solution  [5,  Ch.  4].  One  of  these  categories  is  the  low-energy-coherence 
(LEC)  processes,  in  such  cases  the  energy  is  distributed  over  large  number  of 
coordinates  and  for  which  all  of  the  eigenvalues  are  small  compared  to  the  white 
-  noise  one-sided  power  spectral  density.  It  may  appear  that  the  LEC  condition 
implies  poor  performance  and  is  therefore  uninteresting.  This  is  not  true  because 
the  receiver  output  is  obtained  by  combining  a  large  number  of  components,  and 
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Eqs.  (1.11)  and  (1.12)  are  highly  dependent  on  the  scattering  function  of  the 
channel  through  their  dependence  on  d,  therefore;  in  order  to  implement  the 
optimum  receiver  and  to  evaluate  its  performance  one  needs  to  know  the  scatter¬ 
ing  function, cr(r,f).  Lack  of  complete  knowledge  of  cr(r,f)  will  result  in  degrada¬ 
tion  of  the  receiver’s  performance. 


1.4  PERFORMANCE  OF  A  MISMATCHED  RECEIVER 

Under  mismatch  conditions,  the  square  root  of  the  numerator  in  Eq.  (1.13) 
can  be  written  in  terms  of  the  true  channel  scattering  function,  at(r,f),  and  the 
assumed  known  to  the  receiver  scattering  function,  <ra(r,f),  as  follows 


e//|  H.Wf  I  H0|=^T/7//<rt(r,fK(r'  S  )  |  x(r-r'  }-<  )  |  2drdfdr'  df'(1.15) 


and  the  denominator  is 


Var/ / 1  H.W jL J  /  /K(r,<X(r'  S  )  |  x(r-r  S-t  )  |  2drdfdr'  df'  ,  (1.16) 

V  J  N0  _oo 


where  xih^)  is  the  ambiguity  function  of  the  complex  envelope  of  the  transmitted 
signal,  namely 


00  T  1 

x(t, v)=  f  u(t)u*(t-r)exp  Ij27ri4  dt. 


(1.17) 


-00 


Case  1:  Scattering  Functions  With  Unknown  Parameters 

The  first  case  under  consideration  is  the  case  when  there  is  a  mismatch  in 
the  mean  time  delay  and  the  mean  Doppler  shift  by  say  X!  and  x2  respectively. 
Then  the  true  and  the  assumed  scattering  functions  can  be  related  as  follows 

(1.18) 


^a(r>f)=^t(r-xHf'x2) 


a 


a 


i 

3 
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Substituting  in  Eqs.  (1.15)  and  (1.16),  we  obtain 


/ /  f  )  I  X(r-r  +x„f-f/  +x2)  |  2drdfdr'  df' 

-OO 


1 1  JJvaWMt  j  )  I  x(r-r'  .f-f'  )  |  2drdfdr  df' 

-OO 


(1.19) 


In  order  to  proceed  further  with  the  calculations,  one  needs  to  specify  the 
shape  of  the  true  scattering  function  and  the  envelope  of  the  transmitted 
waveform.  If  the  true  scattering  function  of  the  channel  is  assumed  to  have  the 
following  shape 

<Tt(r’f)=2^BeXP  [-ri/2L2-f2/2B2  ]•  (1-20) 


where  Er  is  the  total  received  energy  when  the  complex  envelope  of  the  transmit 
ted  waveform  has  unit  energy.  And  the  envelope  of  the  transmitted  signal  is 

/  1/x/T  0<t<T 
u(t)  elsewhere. 

Then  the  true  autocorrelation  function  is 


R(t,r)= 


Er/T  exp 
0 


[-2(7rB(r-t))2  j|Q(max(r-T)/L)-Q(min(t,r)/L)|  |  r- 1  |  <T 

elsewhere.  (1*^2) 


Equation  (1.19)  can  be  evaluated  numerically.  Then  the  probabilities  of  detec¬ 
tion,  PD,  and  false  alarm,  PF,  can  be  determined  by  substituting  in  Equations 
(1.11)  and  (1.12).  Figures  1.2  and  1.3  demonstrate  the  change  in  the  probability 
of  detection,  PD,  versus  SNR  for  probability  of  false  alarm,  Pp,  of  10-4  and 
different  values  of  mean  time  delay,  x1?  and  the  mean  frequency  shift,  x2.  It  is 
clear  that  the  receiver  is  more  sensitive  to  misalignment  in  frequency  shifts  than 
misalignments  in  time  delays.  Also,  the  amount  of  degradation  in  performance  for 
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under-spread  channels  (s=0.1)  for  x2>B  and  Xj>5L  is  severe  and  the  system 
designer  should  provide  the  adequate  circuitry  to  insure  synchronization. 

.  Case  2:  Completely  Unknown  Scattering  Functions 

If  the  shape  of  the  true  scattering  function  is  not  known,  the  receiver  has, 
either  to  estimate  it  or  use  an  approximation  for  <rt(r,f).  In  both  cases  the  true 
and  the  assumed  scattering  functions  may  differ  in  shape  and  the  resulting  meas¬ 
ure  of  performance,  d2,  is 


1 1  11° t(r»fWr  f  )  I  X(r-r'  .f-f'  )  |  2  drdfdr  df' 

1  -0° 

dJ=  1  - 

/  /  /  f‘r.(rJMr  /  )  I  x(r~r  ,{-(  )  I  2  drdfdr'  df1 


For  the  present  discussion,  <ra(r,f)  will  be  taken  as 

n_/Er/LB  >f|  -B/2’ 
aa(r»  )  elsewhere, 


(1.23) 


(1.24) 


while  the  true  scattering  function,  tft(r,f)  will  be  represented  by  Eq.  (1.20). 

Figure  1.4  illustrates  the  effect  of  mismatch  in  the  scattering  functions  for  this 
specific  case. 

For  over-spread  channels,  (s  >  10.0),  equation  (1.22)  is  very  closely 
approximated  by  a  knife-edge  (ridge)  function  along  the  t-r  axis.  In  this  case  we 
expect  that  the  system  become  more  sensitive  for  any  time  delays,  frequency 
shifts,  or  mismatch  in  the  shape  of  the  scattering  functions. 

It  is  noticed  also  that  higher  values  of  Pp  will  not  change  the  previous 


results  significantly. 
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1.5  CONCLUSIONS 

In  order  to  implement  the  optimum  receiver  over  fading  dispersive  chan¬ 
nels  in  practice,  one  needs  the  filter  impulse  response  h(t,r).  This  is  obtained  by 
-  solving  the  integral  equation  (8)  which  depends  on  the  scattering  function  o(r,f) 
through  the  autocorrelation  function  R(t,r).  Both  measurement  of  the  <r(r,f)  and 
the  solution  of  the  integral  equation  present  barriers  to  the  actual  receiver  imple¬ 
mentation. 

This  triggered  the  need  to  study  the  sensitivity  of  a  mismatched  receiver. 
We  began  by  studying  the  effects  of  mismatch  in  the  location  of  the  scattering 
function.  It  has  been  demonstrated  that  for  a  Gaussian  shaped  scattering  func¬ 
tions  and  under-spread  channels,  the  amount  of  degradation  is  severe  if  Xj>5L 
and  x2>B.  We  next  studied  the  degradation  due  to  mismatch  in  the  shape  of  the 
scattering  functions.  A  degradation  of  about  1  dB  is  encountered  if  a  rectangular 
scattering  function  instead  of  a  Gaussian  shaped  scattering  function  is  used  at 
the  receiver,  provided  that  the  true  scattering  function  was  Gaussian. 

It  worthwhile  mentioning  that  if  we  change  the  role  of  at  and  the 
amount  of  the  degradation  will  not  the  same.  As  a  matter  of  fact,  it  can  be 
shown  that  if  <ra( r,f)  is  Gaussian  but  <7t(r,f)  is  rectangular,  then  the  degradation  in 
performance  is  less  severe. 

For  over-spread  channels,  the  autocorrelation  functions,  R(t,r)  tend  to  con¬ 
centrate  around  the  t=T  axis  and  any  misalignment  in  location  or  mismatch  in 
shape  will  have  a  greater  degradation  effects  on  the  performance  of  the  receiver. 

The  study  emphasis  the  fact  the  synchronization  over  fading  dispersive 
channels  is  important  and  in  some  cases  the  degradation  due  to  missynchroniza- 
tion  may  be  more  severe  than  mismatch  in  the  shape  of  the  scattering  function, 
this  mismatch  may  be  a  result  of  error  in  measuring  of  the  scattering  function. 
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Figure  1.2  Mismatch  in  the  Location  of  the  Scattering  Function 


Channel  Spread  S=0.10 


SNR  in  dB 

Figure  1.3  Mismatch  in  the  Location  of  the  Scattering  Function 
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CHAPTER  H  | 

SERIAL  SEARCH  SYNCHRONIZATION  OVER  | 

FADING  DISPERSIVE  CHANNELS  1 


Although;  the  engineering  importance  of  fading  dispersive  channels  has  increased 
markedly  in  the  recent  years.  The  problem  of  synchronization  over  such  channels 
has  not  been  given  the  attention  it  deserves.  The  following  analysis  derives  and 
evaluates  upper  and  lower  bounds  on  the  mean  and  variance  of  synchronization 
time  for  a  very  general  serial  search  system.  The  objective  is  to  find  a  simple 
result  that  the  system  designer  will  be  able  to  use  to  make  design  trade-off  stu¬ 
dies  to  minimize  mean  synchronization  time. 

It  has  been  found  that  the  results  are  highly  dependent  on  the  scattering  function 
of  the  channel,  the  degree  of  the  channel  spread,  number  of  the  cells  used  in  the 
search  procedure,  and  the  the  time  required  to  reject  an  incorrect  cell  when  a 
false  alarm  occurs. 

2.1  INTRODUCTION 

Synchronization  is  a  fundamental  problem  in  digital  communication, 
radar,  sonar,  and  navigation  systems.  Power-efficient  receivers  generally  require 
the  existence  of  a  clock  that  is  accurately  time-aligned  with  the  received  pulses, 
and  a  local  carrier  reference,  that  agrees  closely  in  frequency  and  phase  with  the 
received  carrier.  In  general  terms,  two  sequences  of  events  are  said  to  be  synchro¬ 
nous  if  the  corresponding  events  in  the  two  sequences  occur  simultaneously,  with 
one  of  the  two  sequences  of  events  takes  place  at  the  transmitter,  and  the  other 
takes  place  at  the  receiver.  Due  to  transmitter  oscillator  instability  and  propaga¬ 
tion  effects,  the  two  sequences  of  events  may  be  misaligned. 
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Losing  synchronization  reduces  the  efficiency  in  the  data  detection  process; 
because  the  inaccurate  symbol  sync  directly  reduces  the  probability  of  making 
correct  decisions.  Moreover,  when  a  loss  of  synchronization  occurs,  it  *nay  some¬ 
times  lead  to  successive  errors  before  it  is  regained;  these  successive  errors  affect 
the  overall  performance  of  the  system. 

Synchronization  is  defined  simply  as  the  process  of  bringing  about,  and 
retaining,  a  synchronous  situation.  It  is  generally  convenient  to  separate  the  syn¬ 
chronization  process  into  two  distinct  modes.  In  the  first  mode,  the  clock  syn¬ 
chronization  mode,  the  clocks  which  regulate  the  two  sequences  being  synchron¬ 
ized  (i.e.,  the  transmitter  and  receiver  clocks)  are  forced  to  run  at  the  same  rate. 
In  the  second  mode,  the  higher  order  synchronization  mode,  a  corresponding  pair 
of  events  in  the  two  sequences  are  identified  and  made  to  occur  simultaneously, 
and  if  the  sequences  are  progressing  at  the  same  rate,  the  sequences  are,  and  will 
remain,  synchronized. 

If  the  transmitter  and  the  receiver  clocks  are  both  sufficiently  stable  rela¬ 
tive  to  the  required  synchronization  accuracy,  the  clock  synchronization  mode 
may  be  by-passed.  Most  generally,  however,  this  will  not  be  the  case  and  some 
technique  must  be  devised  to  provide  the  needed  clock  synchronization. 

Once  the  transmitter  and  receiver  clocks  have  been  synchronized  the 
second  mode  of  the  synchronization  process  begins,  this  includes  symbol,  code 
word,  data  word,  and  frame  synchronization.  The  second  mode  can  be  further 
subdivided  into  two  components,  the  first  component  is  the  determination  of  the 
initial  parameters  from  whatever  a  priori  information  available.  This  component 
is  known  as  the  acquisition  phase  The  second  component  is  to  maintain  syn¬ 
chronization  after  initial  acquisition  .  This  problem  is  known  as  the  tracking 
phase. 


Although;  the  problem  of  synchronization  over  additive  white  Gaussian 
noise  (AWGN)  channels  has  been  investigated  extensively  [7-10],  synchronization 
over  fading  dispersive  channels  has  not  been  given  the  attention  it  deserves,  in 
.  spite  of  the  fact  that  the  the  engineering  importance  of  such  channels  has 
increased  markedly  in  the  recent  years.  Fading  dispersive  channels  are  usually 
best  described  as  random  linear  time-varying  filters  .  In  many  applications  it  is 
reasonable  to  suppose  that  the  impulse  response  of  the  filter  is  a  sample  function 
of  a  Gaussian  random  process. 

In  most  of  the  work  done  concerning  detection  over  fading  dispersive  chan¬ 
nels,  the  assumption  that  synchronization  is  available  is  either  understood  impli¬ 
citly  or  declared  explicitly. 

2.2  SYNCHRONIZER  STRUCTURE 

In  communicating  over  fading  dispersive  channels,  the  received  energy  is 
peaked  at  some  point  P  (see  Figure.  2.1),  which  is  shifted  from  the  transmitter 
clock  pulse  and  oscillator  frequency,  point  T,  by  7  sec.  in  time  and  f  Hz  in  fre¬ 
quency,  where  7  and  f  are  the  mean  time  delay  and  mean  frequency  shift  of  the 
channel  respectively.  All  the  points  are  referred  to  the  receiver  clock  pulse  and 
oscillator  frequency  which  are  located  at  the  origin  of  the  figure.  It  is  clear  that 
no  matter  what  offset  exists  between  the  receiver  R  and  the  transmitter  T,  the 
important  information  to  the  receiver  is  the  vector  &  ,  where  &  =  (xj,  x2).  There¬ 
fore  the  main  objective  of  the  receiver  is  to  come  up  with  a  good  estimate  of  the 
vector  &• 

The  vector  £  is  the  misalignment  in  time  and  frequency  between  the 
receiver  and  transmitter  is  the  result  of  two  effects,  the  first  is  the  contribution  of 
the  channel,  and  the  second  is  due  to  the  lack  of  synchronization  between  the 


receiver’s  and  transmitter’s  clock  pulses  and  oscillator  frequencies.  This  lack  of 
synchronization  always  exists  initially  even  if  the  channel  introduces  no  delay  in 
time  or  shift  in  frequency.  Since  these  two  effects  are  of  the  same  nature  they  can 
.  be  lumped  together  into  the  mean  values  of  the  scattering  function,  therefore 
from  now  and  on  we  will  assume  that  the  scattering  function  is  of  the  form 

^(r,f:x)=^(r-x1,f-x2).  (2.1) 

For  such  scattering  function  the  correlation  function  of  the  complex  envelope  of 
the  output  of  the  channel  for  single  pulse  is 

R(t,nx)=R0(t-x1r-x1)exp  |j27rx2(r-t)  j  (2.2) 

where  Rc(t,r)  is  the  complex  correlation  function  of  the  complex  envelope  of  the 
output  of  the  channel  in  the  case  of  perfect  synchronization. 

The  theory  of  maximum  likelihood  (ML)  estimation  provides  us  with  a 
technique  to  estimate  the  channel  parameters  needed  for  the  sync  procedure.  In 
[4]  we  derived  the  likelihood  function.  The  rule  is  x.  is  more  likely  than  x  iff 

W(x)  >  W(*  ), 

where 

W(x.)= J  J  z*(t)h(t,r;x.)z(r)dtdr.  (2.3) 

h(t,r:x)  is  the  solution  of  the  following  integral  equation 

/R(t,u:x)h(u,nx)du+N0h(t,r:x=R(t,r:2().  (2.4) 

N0/2  is  the  two  sided  power  spectral  density  of  the  AWGN,  z(t)  is  the  complex 
envelope  of  the  received  waveform  (see  Figure  1.1  ). 

Under  low-energy-coherence  condition  [5]  equation  (2.3)  is  reduced  to 
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W(x)= / / z*(t)R(t,T:x)z(r)dtdr.  (2.5) 


The  precise  implementation  of  the  maximum  likelihood  estimator  would  involve  a 
continuum  of  detectors  (or  calculations)  for  each  value  of  x  in  the  region  of 
interest.  Certainly,  in  general  such  estimator  cannot  be  built.  One  realizable 
approximation  to  the  maximum  likelihood  estimator  results  from  specifying  a 
finite  set  of  points,  say  {  x,  ,i=l,...,N},  and  record  z(t)  then  perform  the  parallel 
processing  operation  shown  in  Figure  2.2  in  which  h(t,nx)  is  viewed  as  a  linear 
time-varying  filter  with  as  the  input  time  variable  and  t  as  the  output  time  vari¬ 
able.  The  output  of  this  processor  is  N  numbers,  each  is  an  output  of  a  detector 
matched  to  a  particular  x*  The  ML  estimator  finds  the  largest  detector  output 
and  assumes  that  the  correct  value  of  x  is  the  one  corresponding  to  that  detector. 
The  result  is  a  parallel  processor,  in  each  parallel  branch  the  likelihood  function 
is  evaluated  and  the  maximum  is  chosen  to  represent  the  maximum  likelihood 
estimate  of  the  vector  of  parameters  under  consideration.  This  synchronizer 
which  is  optimum  in  the  sense  that  it  achieves  synchronization  with  a  given  pro¬ 
bability  in  the  minimum  possible  time  requires  a  detector  for  every  cell  and  thus 
is  not  optimum  in  a  minimum  hardware  sense.  Thus  minimum  acquisition  time 
systems  is  never  implemented  because  of  excessive  hardware  complexity.  A  syn¬ 
chronization  system  that  evaluates  the  cells  serially  until  the  correct  cell  is  found 
is  said  to  use  serial  search.  Such  systems  are  designed  to  achieve  a  compromise 
between  acquisition  time  and  reasonable  complexity  without  compromising  any 
other  important  system  characteristic.  The  performance  of  such  serial  search  syn¬ 
chronizers  are  analyzed  in  the  next  section. 
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2.3  PERFORMANCE  ANALYSIS 

If  the  uncertainty  region  is  divided  into  N  cells  as  shown  in  Figure  2.3, 
where  the  size  of  each  cell  is  proportional  to  the  width  of  the  autocorrelation 
-  ambiguity  function.  This  function  is  the  expected  value  of  the  random  variable 
W(x).  A  complete  discussion  of  the  properties  of  this  autocorrelation  ambiguity 
function  is  given  in  [ll]. 

Because  it  is  equally  likely  that  the  correct  vector  is  in  any  cell,  the  search 
can  begin  at  any  corner  of  the  uncertainty  region.  The  search  will  advance 
through  one  cell  at  a  time  until  N  cells  have  been  evaluated.  If  synchronization 
has  not  been  achieved  at  that  time,  a  retrace  will  start  the  search  over  again  at 
I  the  starting  position. 

\  The  mean  synchronization  time  is  calculated  by  considering  all  possible 
sequences  of  events  leading  to  a  correct  synchronization.  An  event  in  the  proba¬ 
bility  space  being  considered  is  defined  by  a  particular  location,  n,  for  the  correct 
cell,  a  particular  number  of  missed  detection,  j,  of  the  correct  cell,  and  a  particu¬ 
lar  number  of  false  alarms,  k,  in  all  incorrect  cells  evaluated.  The  total  synchroni¬ 
zation  time  for  a  particular  event  defined  by  (n  j,k)  is 

T(n  j,k)  =  nTe  +  jNTe  +  kTfa,  (2.6) 

where  Te  is  the  evaluation  time  for  each  cell,  and  Tfa  is  the  time  required  to 
reject  an  incorrect  cell  when  a  false  alarm  occurs.  It  can  be  shown  that  the  pro¬ 
bability  of  the  event  (n,j,k)  is 

pr(n,j,k)  =  -Jj-pd  (1-Pd y  (£ )p^(l-pfa)K_k.  (2.7) 

The  mean  synchronization  time  is 

T8=  £T(n,j,k)pr(n,j,k).  (2.8) 

n,j,k 
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After  some  straightforward  algebraic  manipulation  [30-Ch.  10],  equation  (12)  is 
reduced  to 

T,  =  (N-l)Trf(^?i)+^,  (2.9) 

2Pd  Pd 

where 

Tad  ~  Te  +  TfaPfa,  (2.10) 

is  the  average  dwell  time  on  each  cell.  Equation  (2.9)  also,  has  been  derived  in 
[12]  using  signal  flow  graph  techniques.  It  can  be  shown  also,  that  the  variance 
of  the  synchronization  time  is  approximated  by 

«  T^N2(i--i-+  J_)  (2.11) 

12  Pd  Pd 

The  approximation  is  valid  for  N>>1,  pfa<<l,  and  l-pd<<l. 

Calculations  of  pd  and  pJa: 

The  serial  search  synchronizer  under  consideration,  evaluates  each  cell  by 
estimating  whether  or  not  signal  energy  is  present  at  the  output  of  each  processor 
in  Figure  2.2  or  simply 

W(x,)>7.  (2.12) 

H0 

where,  Hx  is  the  hypothesis  that  the  x,  is  the  correct  while  H0  is  the 
hypothesis  that  Xj  is  not  the  correct  cell.  This  problem  is  identical  to  the  prob¬ 
lem  of  detection  over  fading  dispersive  channels  which  is  well  treated  in  |3]  and 
[5].  Van  Trees  [5,  Ch.  4]  showed  that 


where 


and 


00 

Q(*)=f-^exP 


(2.15) 


(2.16) 


In  order  to  proceed  further  with  the  calculations,  one  needs  to  know  the 
shape  of  both  the  channel  scattering  function  and  the  complex  envelope  of  the 
transmitted  wave  form.  Assuming  Gaussian-shaped  functions 

o(r,f:x)  =  exp  [-(r-x1)2/2L2-(f-x2)2/2B2  ],  (2.17) 


where  Er  is  the  total  received  energy  when  the  complex  envelope  of  the  transmit¬ 
ted  waveform  has  a  unit  energy.  The  complex  envelope  of  the  transmitted 
waveform  is 

u(t)  =  (2— )1//4exp  [-at2  1.  (2.18) 

7 r  t  J 


It  has  been  shown  that  [4]  the  optimum  value  of  a  is 


(2.19) 


Figures  2.4  and  2.5  demonstrate  how  pd  behave  as  a  function  of  the 
signal-to-noise  ratio  (SNR)  and  pfa  for  both  under-spread,  and  over-spread  chan¬ 
nels  respectively.  These  results  are  valid  only  under  LEC  condition,  and  if  the 
width  of  each  cell  is  greater  than  or  equal  to  the  width  of  the  autocorrelation 
ambiguity  function  in  both  time  and  frequency.  This  put  a  restriction  on  the 
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accuracy  of  the  estimate,  and  the  number  of  cells  N 

Tfl 


n^ap 


(24) 


where  A  and  T  are  the  width  of  the  autocorrelation  in  time  and  frequency  respec¬ 
tively. 


2.4  CONCLUSIONS 

A  method  to  estimate  the  mean  synchronization  time  for  synchronization 
over  fading  dispersive  channels  has  been  shown.  Equation  (13)  demonstrates  the 
dependence  of  the  mean  synchronization  time,  Tg  on  pd,  p(a,  Te,  T{a,  and  N.  To 
minimize  Ts  one  has  to  select  the  optimum  set  of  parameters.  Some  of  the  these 
parameters,  as  N,  can  be  controlled  by  the  designer,  the  remaining,  pd,  pfa,  and 
Te  are  dependent  on  the  shape  of  the  scattering  function,  the  amount  of  spread 
of  the  channel,  and  SNR. 
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CHAPTER  ffl 


CLOSED  LOOP  SYNCHRONIZER  OVER 
FADING  DISPERSIVE  CHANNELS 

In  this  chapter  an  optimum  closed  loop  structure  for  symbol  synchronizer 
over  fading  dispersive  channels  is  developed.  The  optimum  synchronizer  is  simi¬ 
lar  to  those  already  known  symbol  synchronizers  that  are  being  used  over  addi¬ 
tive  white  Gaussian  noise  (AWGN)  channels,  except  that  the  control  signal  is 
completely  random;  even  in  the  absence  of  the  AWGN.  The  reason  is  clear,  and 
simply  because  of  the  structure  of  the  received  signal  over  fading  dispersive  chan¬ 
nels.  It  is  shown  that  the  synchronizer  structure  ‘is  highly  dependent  on  the 
scattering  function  of  the  channel.  For  a  Gaussian  shaped  scattering  function,  a 
simple  closed  loop  structure  is  obtained  and  expressions  for  tracking  error  statis¬ 
tics  are  derived.  The  dependency  of  such  statistics  on  the  signal-to-noise  ratio 
and  on  the  spread  of  the  channel  is  investigated. 

3.1  INTRODUCTION 

It  is  generally  convenient  to  separate  the  synchronization  process  into  two 
distinct  modes.  In  the  first  mode,  the  clock  synchronization  mode,  the  clocks 
which  regulate  the  two  sequences  being  synchronized  (i.e.,  the  transmitter  and 
receiver  clocks)  are  forced  to  run  at  the  same  rate.  In  the  second  mode,  the 
higher  order  synchronization  mode,  a  corresponding  pair  of  events  in  the  two 
sequences  are  identified  and  made  to  occur  simultaneously,  and  if  the  sequences 
are  progressing  at  the  same  rate,  the  sequences  are,  and  will  remain,  synchron¬ 
ized. 

If  the  transmitter  and  the  receiver  clocks  are  both  sufficiently  stable 
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relative  to  the  required  synchronization  accuracy,  the  clock  synchronization  mode 
may  be  by-passed.  Most  generally,  however,  this  will  not  be  the  case  and  some 
technique  must  be  devised  to  provide  the  needed  clock  to  maintain  synchroniza- 
_  tion  after  initial  acquisition  .  This  problem  is  known  as  the  tracking  phase. 

The  advent  of  highly  compact,  inexpensive  digital  computers  has  now  made  it 
possible  to  exploit  well-known  results  from  statistical  estimation  theory  and  con¬ 
trol  theory  to  develop  systems  that  automatically  respond  to  changing  signal 
environment.  This  self  adjusting  or  adaptive  capability  renders  the  operation  of 
such  systems  more  flexible,  reliable,  and  more  importantly  offers  improved  recep¬ 
tion  performance  that  would  be  difficult  to  achieve  in  any  other  way.  These  rea¬ 
sons  motivate  the  idea  of  developing  a  synchronizer  which  updates  itself  (closed 
loop  structure). 

In  the  following  sections,  a  closed  loop  structure  synchronizer  is  derived 
and  its  performance  is  investigated. 

As  in  chapter  II,  the  vector  x*.  is  used  to  indicate  the  misalignment  in  time 
and  frequency  between  the  receiver  and  transmitter  is  the  result  of  two  effects, 
the  first  is  the  contribution  of  the  channel,  and  the  second  is  due  to  the  lack  of 
synchronization  between  the  receiver’s  and  transmitter’s  clock  pulses  and  oscilla¬ 
tor  frequencies.  This  lack  of  synchronization  always  exists  initially  even  if  the 
channel  introduces  no  delay  in  time  or  shift  in  frequency.  Since  these  two  effects 
are  of  the  same  nature  they  can  be  lumped  together  into  the  mean  values  of  the 
scattering  function.  This  will  result  in  scattering  functions  as  described  by  equa¬ 
tion  (2.1).  The  optimum  estimate  of  x  is  such  that 

W(x)  >  W(x'  )  (3.1) 

where 


IWUB  IJglUMW  I  JP  W  W  iff  Hf Aw  g' M  W  KMWW  **  *-" HgF  VmsftwwW.wm'^WlWi’JlWlJWi 


48 

b(t,r:x)  is  the  solution  of  the  following  integral  equation  given  by  equation  (2.4), 
and  Nc/2  is  the  two  sided  power  spectral  density  of  the  AWGN,  and  z(t)  is  the 
complex  envelope  of  the  received  waveform  (see  Figure  1.1).  Under  low-energy- 
coherence  condition  [5]  equation  (3.2)  is  reduced  to 

w(x)=-^j//z‘(t)R(t,T:2c)z(r)dtdr.  (3.3) 

The  low-energy-coherence  condition  is  more  realistic  for  over-spread  channels. 

In  the  case  of  repeated  transmission  the  transmitted  waveform  can  be 
written  as 

s(t)=X)s1(t-(i-l)/p)  (3.4) 

1=1 

where  p  is  the  transmission  rate.  It  can  be  shown  that  [4]  the  likelihood  function 
for  estimating  Xt  in  the  case  of  repeated  transmission  is 

W(x)=-jy  £  //zi*(t)R|(t,nx.)z,(r)dtdr,  (3.5) 

where  R,(t,7:x)  is  the  complex  correlation  function  associated  with  the  i-th 
transmitted  waveform. 

As  mentioned  before,  the  precise  implementation  of  the  maximum  likeli¬ 
hood  estimator  would  involve  a  continuum  of  detectors  (or  calculations)  for  each 
value  of  x.  in  the  region  of  interest.  Certainly,  in  general  such  estimator  cannot  be 
built.  One  realizable  approximation  to  the  maximum  likelihood  estimator  results 
from  specifying  a  finite  set  of  points,  say  {  &  ,i=l,...,N},  and  record  z(t)  then 
perform  the  parallel  processing  operation  in  which  h(t,r.x)  is  viewed  as  a  linear 
time-varying  filter  with  r  as  the  input  time  variable  and  t  as  the  output  time 
variable.  The  output  of  this  processor  is  N  numbers,  each  is  an  output  of  a  detec¬ 
tor  matched  to  a  particular  x*  The  ML  estimator  finds  the  largest  detector  output 
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and  assumes  that  the  maximum  likelihood  estimate  of  2S*  is  the  one  corresponding 
to  that  detector.  The  result  is  a  parallel  processor,  in  each  parallel  branch  the 
likelihood  function  is  evaluated  and  the  maximum  is  chosen  to  represent  the 
-  maximum  likelihood  estimate  of  the  vector  of  parameters  under  consideration. 
This  synchronizer  which  is  optimum  in  the  sense  that  it  achieves  synchronization 
with  a  given  probability  in  the  minimum  possible  time  requires  a  detector  for 
every  cell  and  thus  is  not  optimum  in  a  minimum  hardware  sense.  Thus 
minimum  acquisition  time  systems  is  never  implemented  because  of  excessive 
hardware  complexity.  A  synchronization  system  that  evaluates  the  cells  serially 
until  the  correct  cell  is  found  is  said  to  use  serial  search.  Such  systems  are 
designed  to  achieve  a  compromise  between  acquisition  time  and  reasonable  com¬ 
plexity.  The  performance  of  a  serial  synchronizer  is  discussed  in  [13]. 

In  the  following  sections  a  closed  loop  synchronizer  is  derived  and  its  per¬ 
formance  is  analyzed. 

3.2  CLOSED  LOOP  STRUCTURE 

Our  objective  is  a  closed  loop  structure  that  can  be  derived  from  the  ML 
equations.  The  reluctance  to  develop  a  closed  loop  structure  for  the  synchronizer 
is  due  principally  to  the  continuously  changing  nature  of  the  characteristics  of 
the  channel.  The  technique  we  will  use  to  transform  the  open  loop  structure 
(parallel  processor)  derived  from  the  maximum  likelihood  theory  to  the  desired 
closed  loop  structure  is  the  gradient  approach.  The  gradient  approach  is  very 
popular  since  it  is  a  relatively  simple  and  generally  well  understood  method  that 
permits  the  solution  of  a  large  class  of  problems.  When  the  likelihood  function  is 
near  quadratic,  then  the  performance  measure  can  be  visualized  as  a  bowl¬ 
shaped  surface,  so  the  synchronizer  has  the  task  of  continually  seeking  the  “top 


of  the  bowl”.  In  this  case,  seeking  the  maximum  of  the  likelihood  function  can  be 
accomplished  by  “hill  climbing”  method  of  which  the  various  gradient  methods 
are  representative.  For  the  ML  function  (3.5),  the  gradient  is  obtained  by 
differentiating  with  respect  to  the  vector  i  to  yield 

j_,,2  (3.6) 

UXl  i=l  OTj 

In  order  to  proceed  further  with  the  analysis,  one  needs  to  know  the  shape  of  the 
channel  scattering  function.  By  assuming  a  Gaussian- shaped  scattering  function 
of  the  form 


^r»f:20=  2^b  exp  [^r"xi)2/2L2-(f-x2)2/2B2  ], 


where  Er  is  the  total  received  energy  when  the  complex  envelope  of  the  transmit¬ 
ted  waveform  has  unit  energy.  Using  such  scattering  function,  equations  (12)  can 
be  reduced  to 


^X|W(2.)-  £  //z,,(t)Rl,(t,ri)2,('-)dtdT,  j— 1,2 


(3.8)) 


where 


Rii(t,7:20= / / ro(r,f :&)u,(t-r )u, *(r-r)exp  |j27rf(r— t)  jdrdf,  (3.9-a) 

and 

Ri2(t,r:&)= J  J  fa^fiiju^t-rju^r-rjexp  p27rf(r-t)  jdrdf.  (3.9-b) 

Equating  (3.0-a)  and  (3.0-b)  to  zero  and  solving  the  two  transcendental  equations 
simultaneously  will  yield  the  maximum  likelihood  estimate  of  &. 


An  alternative  to  be  considered  in  this  paper  is  to  use  the  partial  derivatives 


^  J-1'2 


(3.10) 
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to  control  the  current  estimate  of  Since  this  expression  will  be  zero  when  & 
equals  the  true  vector  and  since  it  is  monotonically  increasing  function  of  x  in 
the  vicinity  of  2L=Xt>  a  feedback  device  should  be  able  to  force  &  to  converge  to 
-  the  maximum  likelihood  estimate  of  Xt-  The  resulting  symbol-tracking  device  is 
shown  in  Figure  3-1.  The  output  of  the  box  labeled  “accumulator”  is 

ej  =  XjW(x)  -  Wj(x),  j=l,2;  (3.11) 

where 

Wj&)  =  £//zi*(t)Rij(t,7:x)zi(r)dtdn  j=l,2;  (3.12) 

i=l 

The  error  signals  ej,  j=l,2  are  random  variables,  their  means  and  variances  are 
given  by 


E^ej  ^=XjE^  W(x)  Wj(x)  ^ 

Var^e,  |=xJ2var|  W(x)  |+var  j  Wj(x)|-2XjCov|  W,Wj  | 


(3.13) 

(3.14) 


3.3  NUMERICAL  RESULTS 

To  evaluate  the  different  variances  one  needs  to  specify  the  complex 
envelope  of  the  transmitted  signal.  For  Gaussian-shaped  scattering  function,  it 
has  been  found  that  a  Gaussian-shaped  envelope  in  the  form 

u(t)  =  (2a/7r)1/4exp  [-at2  j,  (3.15) 

is  optimum  from  the  detection  point  of  view  [4].  Other  reasons  for  the  choice  of 
(3.7)  and  (3.15)  for  the  forms  of  the  scattering  function  and  the  complex  envelope 
of  the  transmitted  waveform  are  based  on  the  fact  that  the  Gaussian  function 
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1.  is  a  simple  and  smooth  function, 

2.  has  elegant  properties  when  integrated  against  another  Gaussian  function 
on  the  infinite  interval, 

3.  can  be  used  to  approximate  many  of  the  finite-duration  pulse  signals  com¬ 
monly  used  in  communications,  and 

4.  is  easily  generated  [14]. 

It  has  been  shown  that  the  value  of  a  in  (3.15)  which  minimizes  the  area  of 
the  uncertainty  ellipse  [4]  is 


a  =  ir- 


B 


(3.16) 


The  statistics  of  W(&)  and  W:  are 


E 


|w|=E  JsNR+-i-//R,'(t,r:jtl)Ri(t,rjOdtdr, 

E{W,|=  E  k  SNR+^L.//R|,(t,rat)Rij(t,rat)dtdr  ,  (3.18) 

{whl,  ‘(t.rj^jR^t^^Ri^t^nxtJR^tj.^gcJdtdrdtidr! 


f^///R**(t’ri:^t)Ri(t'T:i)Ri(r*ri:i)dtdrc,7'i 


+TJj///Ri*(tl»T:2£t)Ri(tl»t:i)Ri(t*r:i)dtdrdtl 
+  T77  //Ri*(t»T^)Ri(t’T:2L)dtdr  j  (3.19) 

the  quantity  var/wj\  has  a  similar  expression  as  in  equation  (3.19) 
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expect  Rj(t ,rx)  is  replaced  by  Rjj(t,r:&) 

cov|w,Wj  |=5]  |^7////Ri*(t,r1:2£t)Ri(t,nx)Ri*(t1,r;xt)Rij(t1,r1:xt)dtdrdt1dr1 

+'^///Ri*(t’rl:it)Ri(t>r:^)Rii(7'»rl:X)dtdTdrl 

+7T?///Ri*(ti>T:Xt)Rij(ti't:20Ri(t>r:a0dtd7'dfci 

f /Rj*(t,r;x)Rjj(t,r:x)dtdr j  (3.20) 

Substituting  in  Equations  (3.13)  and  (3.14),  we  obtain  Figures  3.2,  3.3,  3.4. 
Without  loss  of  generality  we  will  assume  that  x2t=0.  Figure  3.2  illustrates 
the  behavior  of  the  expected  value  of  the  error  signal  as  a  function  of 
6=(x!-xlt)/L  for  two  values  of  SNR  and  channel  spread  of  10.  Observe  that 
there  is  a  range  of  b  for  which  the  expected  value  of  the  error  signal  is 
almost  linearly  related  to  6.  This  region  may  be  selected  as  a  normal 
operating  region  for  the  tracking  loop. 

Figure  3.3  demonstrate  the  behavior  of  the  expected  value  of  the  error  sig¬ 
nal  for  channels  with  spread  =1.0.  Figure  3.4  shows  the  how  the  variance 
of  the  control  signal  changes  with  6  for  different  values  of  SNR. 

3.4  CONCLUSIONS 

A  closed  loop  structure  for  synchronization  over  fading  dispersive  chan¬ 
nels,  with  Gaussian  shaped  scattering  function,  is  developed.  The  effects  of 
various  system  parameters  (such  as  signal-to-noise  ratio,  the  spread  of  the 
channel  and  the  duration  of  the  envelope  of  the  transmitted  waveform)  on 
the  shape  of  the  S-curve  are  demonstrated.  Dependency  of  the  variance  of 


the  tracking  error  on  the  SNR  and  the  spread  of  the  channel  is  investi¬ 
gated. 

It  is  known  that  in  AWGN  environment  the  optimum  tracking  discrimina¬ 
tor  for  arbitrary  wideband  signal  is  a  multiplier  which  forms  the  product  of 
the  received  signal  plus  noise  and  the  first  derivative,  with  respect  to  the 
parameter  under  consideration,  of  the  receiver  generated  replica  of  the 
transmitted  signal.  This  discriminator  is  optimum  in  that  its  output  is  the 
maximum  likelihood  estimate  of  the  that  parameter  difference  between  the 
two  wideband  signals.  In  synchronization  over  fading  dispersive  channel  it 
is  shown  that  the  optimum  tracking  synchronizer  is  a  multiplier  which 
forms  the  product  of  the  received  signal  and  a  filtered  version  of  it,  the 
filter  impulse  response  is  the  first  derivative  with  respect  to  time  (or  fre¬ 
quency)  of  the  autocorrelation  function  of  the  signal  part  of  the  received 
waveform. 


Ri(k.™)  - *(x) - /(.) dt 


Figure  3.1  Closed  Loop  Synchronizer  Structure 


Channel  Spread  S=10.0 


Figure  3.2  Expected  Value  of  the  Error  Signal 


Channel  Spread  S=1. 


Figure  3.3  Expected  Value  of  the  Error  Signal 
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Figure  3.4  Variance  of  the  Error  Signai 


CHAPTER  IV 


SYNCHRONIZATION  OVER  FADING  DISPERSIVE  CHANNELS 
USING  STOCHASTIC  APPROXIMATION  METHODS 

In  this  chapter  we  consider  the  problem  of  synchronization  over  fading 
dispersive  channels.  Using  stochastic  approximations  methods;  a  recursive 
estimation  procedure  is  developed  to  estimate  the  two  parameters  needed 
by  the  synchronizer.  The  result  is  a  closed  loop  in  which  the  component  of 
the  error  signal  is  proportional  to  the  expected  value  of  the  derivative  of 
the  likelihood  function  with  respect  to  the  appropraite  parameter.  This 
method  possesses  the  simple  computational  structure  of  the  stochastic 
approximation  methods,  and  under  certain  regularity  conditions,  it  can  be 
shown  that  the  variance  of  the  error  in  the  estimates  approaches  the 
Cram4r-Rao  bound.  This  chapter  is  organized  as  follows:  In  the  next  sec¬ 
tion  an  introduction  to  the  problem  is  presented,  in  section  4.2  a  brief  dis¬ 
cussion  to  the  method  of  stochastic  approximation  is  given.  In  section  4.3 
the  method  of  stochastic  approximation  is  applied  to  estimate  recursively 
the  parameters  under  consideration.  Lastly  an  example  is  given  in  section 
4.4  to  illustrate  that  the  suggested  procedure  is  asymptotically  efficient  and 
derive  the  conditions  under  which  asymptotic  efficiency  is 

4.1  INTRODUCTION 

In  this  chapter  the  channel  model  used  is  the  same  one  used  in  section  1.2, 
and  the  main  objective  is  to  develop  a  synchronizer  which  updates  itself 
recursively.  In  [15]  Sakrison  described  an  efficient  recursive  estimation  pro¬ 
cedure  to  estimate  some  target  parameters  from  repeated  observations  using 
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stochastic  approximation  methods.  Under  certain  conditions  these  methods 
provide  many  computational  advantages,  and  it  can  be  shown  that  the 
error  in  the  resulting  sequence  of  the  estimates  approaches  the  Cramer-Rao 
bound. 

In  the  following  section  a  brief  coverage  of  the  method  of  stochastic 
approximation  is  given.  This  method  is  then  applied  to  the  problem  of 
estimating  the  vector  &  described  in  section  2.2. 

4.2  THE  METHOD  OF  STOCHASTIC  APPROXIMATION 
Stochastic  approximation  methods  are  applicable  to  any  problem  that  can 
be  formulated  as  some  form  of  regression  problem  in  which  repeated  obser¬ 
vations  are  made.  To  be  more  specific,  let  the  length  of  the  processing 
subintervals  be  denoted  by  T.  Denote  the  data  observed  on  the  kth  subin¬ 
terval  by  zk.  If  our  observations  start  at  t=0  and  the  observed  quantity  is 
a  continuous-time  random  process,  z(t),  then  zk  represents  the  sample  of  the 
process  of  T  sec  duration,  z(t),  (k-l)T<t<kT.  The  objective  is  to  find  the 
value  of  an  unknown  parameter  x  which  solves  the  vector  valued  equation 

mU)  A  E|l(zk,x)j  =  mo  (4.1) 

Robbins  and  Monro  [16]  originally  studied  an  iterative  procedure  of  solving 
this  problem  from  repeated  independent  observations  in  which  f,  m,  and  z 
are  scalar  valued.  A  multidimensional  version  of  this  procedure  is  described 
in  Sakrison  [17]  where  it  is  shown  that  the  sequence  of  estimates  {x*} 
chosen  according  to  the  recursion  relation 

ik+l=!S»  +  'k[l(ZkJSk)-mo]  (4-2) 


61 


converges  to  the  solution  vector  X*  in  the  mean  square  sense, 

E{llXk-XtH2}  =  k^“E|  x)(xk,rxt,.)2|  =  0  (4.3) 

provided  that  the  following  assumptions  are  satisfied. 

1.  There  exist  constants  ao  and  al5  O<a0<a1<oo,  such  that 

a©llx~Xtll2  <  (i-Xt)T  (m(i)-mo)  <  ailli-itll2-  (4-4) 

where  xT  represent  the  transpose  of  the  vector  £. 

2.  The  random  entities  Zi,i=l,2,...  are  identically  distributed  and  statisti¬ 
cally  independent. 

3.  For  all  values  of  £ 

var||lf(zk,2c||2  j<oo.  (4.5) 

4.  The  sequence  of  constants  {ck},  k=l,2,...  are  positive  monotone  decreas¬ 
ing,  and  satisfy 

00  0° 

Eck=°°»  Eck<°°-  M) 

k=  1  k=l 

The  usage  of  the  stochastic  approximation  methods  yields  the  following  advan¬ 
tages: 

1.  Only  a  small  interval  of  data  needs  to  be  processed  at  a  time. 

•  2.  Only  simple  computations  are  required,  even  when  the  functional  depen¬ 

dence  of  the  regression  function  on  the  parameters  of  interest  'is  nonlinear. 

3.  The  method  may  be  employed  in  the  absence  of  the  detailed  knowledge  of 
the  process  statistics  and  in  the  absences  of  the  detailed  knowledge  of  the 
relationship  between  the  desired  parameters  and  the  observed  data. 
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If  sufficient  a  priori  knowledge  concerning  the  statistics  of  z(t)  and  the  functional 
relationship  between  the  parameters  and  the  observed  data  is  available,  the  third 
advantage  can  be  replaced  by  the  following  desirable  property:  the  methods  can 
-.be  made  to  be  asymptotically  efficient. 

We  can  apply  these  results  to  the  synchronization  problem  under  con¬ 
siderations  mainly  to  find  an  asymptotically  efficient  Robbins-Monro  procedure 
for  estimating  jc  from  the  observations,  zk,  k=l,2,....  The  motivation  for  this  goal 
is  clear;  such  a  method  would  be  computationally  simple  and  yet  would  have  per¬ 
formance  which,  for  a  large  number  of  observations,  k,  could  not  be  surpassed 
by  any  other  method,  no  matter  how  computationally  complex. 

To  see  how  we  might  approach  this  objective,  define 

^k(i)=grad|wk(x)|  (4.7) 

where  Wk(x)  is  W(*)  defined  in  (2.3)  based  on  z(t),  (k-l)T<t<kT.  The  gradient 

in  (4.7)  is  with  respect  to  the  partial  derivatives  ,  j=l,2.  We  notice  that 

*Xj 

3^(x)  is  a  random  vector,  even  in  the  absence  of  the  AWGN,  and  a  useful  charac¬ 
terization  is  the  its  expected  value.  Define  the  vector-valued  function  jn(x)  as 

fflkGO— e|iWx)J-  (4.8) 

We  also  notice  that 


I  =0 

“k&Vo  in  general  <4'9> 

Thus  if  Eq.  (4.9)  has  only  the  unique  solution  we  can  hope  to  carry  out 

Robbins-Monro  method  for  estimating  by  using  the  j^<(x.)  of  Eq.  (4.7)  in  the 
procedure. 


The  most  rapid  convergence  is  obtained  when  the  weighting  constants  ck  are 
chosen  to  be  of  the  form  c/n.  It  can  be  shown  that  [17]  if  c=  G-1^),  where  the 
ijth  entry  of  the  matrix  Gfo)  is  given  by 


(4.10) 


then  the  Robbins-Monro  method  is  asymptotically  efficient.  The  only  problem  is 
how  to  pick  c=G-1(Xt)  ?  since  is  unknown.  One  solution  is  to  make  the  pro¬ 
cedure  adaptive  by  substituting  an  estimate  for  G^).  The  “adaptive”  Robbins- 
Monro  procedure  takes  the  form 


Sk+i=Xk+(l/k)G  KtiMx.)- 


(4.11) 


Sakrison  [18]  considered  the  convergence  of  this  method.  Under  certain  assump¬ 
tions,  Sakrison  has  shown  that  C^,  the  covariance  matrix  of  the  error  (x^-Xt), 
satisfies 


QbC*  b)  <  (l/k)(b,Gb)  +  d/k1+7,  d<oo,  7>0  (4.12) 

where  (2L>x)  is  the  inner  product  of  the  in  the  M-dimensional  Euclidean  space,  and 
b  is  an  arbitrary  vector.  That  is,  the  sequence  of  estimates  generated  by  Eq. 

(4.11)  is  asymptotically  efficient.  The  interested  reader  is  referred  to  [18]  for  the 
proof  of  the  last  statement. 


4.3  NUMERICAL  RESULTS 

Combining  Eqs.  (2.3)  and  (4.10),  yields 
TT 

mk  j==-L.ff  ^’(t.rg^J-N^t-r)  j  R1(t,ra)dtdr  j=l,2;  (4.13) 


KWWW TO! 
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where  superscript  denote  partial  differentiation  with  respect  to  the  corresponding 
component  of  x,  In  order  to  proceed  further  we  need  to  specify  the  shape  of  the 
scattering  function,  cr(r,f).  Let 

=  2^3  exP  [“( r-x l )V 2L2-(f-x2 )2/2B2  ],  (4.14) 

also  assume  that  the  complex  envelope  of  the  transmitted  signal  be 

u(t)  =  (-^-^exp  [-at2  ]  (4.15) 

Substituting  in  Eq.  (1.3)  yields 


-|((t-xi)2  +  (r-x2)2  +  2aL2(t-r)2) 


65 


IT! 


IfUWiPJPUl 


and 


Q  = 


SNR2 

0 


exp  [-T  (A*l2  +  AX22)  ]. 


(4.21) 


The  vector  m(x)  is 


m(x)  =  -2  Q  i 


AX! 

Ax2 


(4.21) 


From  (4.21)  and  (4.19),  the  quadratic  form  q  =  -  (i  -  xJT  G(i)  m(x)  is 
q  =  8tQ2|lAx12('72Ax,2  +  t2Ax22  -  k)  +  BAx22('y2Ax22  +  n^Axj2  -  7)  j  (4.22) 


4.4  CONCLUSIONS 

3W 

The  ML  solution  is  the  one  for  which  — —  =  0  In  the  case  of  fading 

ax, 

dispersive  channels,  the  decision  variable,  even  in  the  absence  of  the  AWGN,  is  a 
random  variable,  a  useful  characterization  is  the  expection  of  W.  One  way  to  do 
approximate  this  solution  is  to  update  the  estimate  by  a  quantity  proportional  to 
the  error.  We  wrote  the  error  function  as 

e(i)=E|e(i)|  +  (e(i)  -  E|e(i)|) 

The  first  term  is  used  to  control  the  new  estimate  and  the  second  term  is  treated 
as  an  internisic  noise.  Eq.  (4.11)  simply  indicates  that  we  update  the  estimate  of 
any  component  of  j  by  using  a  quantity  proportional  to  the  error  .  The  propor¬ 
tionality  factor  is  chosen  to  speed  up  the  convergence  process.  The  resulting  The 
estimate  £.  is  asymptotically  efficient  as  long  as  lies  in  a  given  region,  which 
we  called  the  pull-in  region. 
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